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ABSTRACT 

The desorption of molecular species from ice mantles back into the gas phase in 
molecular clouds results from a variety of very poorly understood processes. We have 
investigated three mechanisms; desorption resulting from H2 formation on grains, di- 
rect cosmic ray heating and cosmic ray induced photodesorption. Whilst qualitative 
differences exist between these processes (essentially deriving from the assumptions 
concerning the species-selectivity of the desorption and the assumed threshold ad- 
sorption energies, E t ) all three processes are found to be potentially very significant 
in dark cloud conditions. It is therefore important that all three mechanisms should 
be considered in studies of molecular clouds in which freeze-out and desorption are 
believed to be important. 

Employing a chemical model of a typical static molecular core and using likely 
estimates for the quantum yields of the three processes we find that desorption by H2 
formation probably dominates over the other two mechanisms. However, the physics 
of the desorption processes and the nature of the dust grains and ice mantles are very 
poorly constrained. We therefore conclude that the best approach is to set empirical 
constraints on the desorption, based on observed molecular depletions - rather than 
try to establish the desorption efficiencies from purely theoretical considerations. Ap- 
plying this method to one such object (L1689B) yields upper limits to the desorption 
efficiencies that are consistent with our understanding of these mechanisms. 

Key words: astrochemistry - molecular processes - stars:formation - ISM: abun- 
dances - ISM: dust - ISM: molecules. 



1 INTRODUCTION 

In cold dark interstellar clouds, heavy molecules accumu- 
late onto dust grains, forming icy mantles on their surfaces. 
This process, known as freeze-out, occurs on a timescale 
of a few times 10 9 n^j 1 years (where nu = w(H) + 2n(H2) 
is the total hydrogen nucleon number density) in the ab- 
sence of desorption. This is much less than the expected 
lifetime of a typical molecular cloud so, if freeze-out was 
unlimited, we would expect the majority of observations 
to show no evidence for heavy gas phase species. However, 
observations of molecules such as CO in dark clouds, for 
example L1689B TMC1-CP and L134 N l|Lee at al l 120031 ; 
ISmith et all 12004 : IWakelam et all 120061 ). indicate that des- 
orption processes must be operating for mantle growth to 
be limited. 

Although this conclusion has been accepted for over 
20 years, it is still not fully understood how this des- 
orption occurs. Many possible mechanisms have been pro- 
posed, most of which require impulsive heating of grains 
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whic h can be caused by (a) direct impact of cosmic 
llHasegawa fc Herbstll 19931: lHartauist fc Williamslll990l : 



rays 1 



iLeger et al.lll985l ). (b)5Travs (|Leger et al.ll 19851 ), (c) ultra- 
violet ph otons induced by cosmic rays (cosmic ray photo des- 
orption) (|Hartquist fe Williams! Il990l : iDulev et al.(ll989l ). or 
(d) exothermic reactions occurring on the grain surface 
jAllen fc Robinson! Il975l ; iGarrod et al]|2007l). in particular 
the formation of molecu lar hydrogen I Willacy et al.lll994al ; 
IDulev fc Williams! 1 1993h . Molecules can either be desorbed 
by classical evaporation (ILeger et al.lll985l) or by chemical 
explosions ( Shen et ail |2004| ; IShalabiea fc Greenbersl Il994l : 
ILeger et al.lll985l ). These chemical explosions can only occur 



if the mantle has previously been irradiated with ultra-violet 
radiation which creates radicals. If the grain temperature 
is raised to ~ 27 K, the radicals become mobile and then 
undergo explosive reactions capable of expelling the entire 
grain mantle. We will not consider chemical explosions in 
this paper because in dark clouds with Av > 5 magnitudes it 
is unlikely that the grain mantles w ill have received sufficient 
U V irradiation (Le ger et al.ll 19851) , and any radicals that do 
form are likely to be hydrogenat ed due to the high abu n- 
dance of hydrogen atoms present l|Willacv fc Millar1ll998h . 
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Recent chemical models of dark clouds tend to in- 
clude only desorption via thermal evaporation (which is 
negligible for dark clouds with temperatures of 10 K) 
and/or direct heating by cosmic r ays, u sing the formula- 
tion given bvlHasegawa fc Herbstl (1 19931) (he reafter HH93) 
l|Ruffle fc Herbstl l2000l : iRoberts et al l 120041 ). This paper, 
therefore, aims to test the assumption that desorption via 
direct cosmic ray heating is the only effective non-thermal 
desorption mechanism relevant to dark molecular cores. In 
this study we adopt the model of a molecular cloud as being 
composed of an ens emble of dense (da rk) cores in a more 
diffuse background (Ga rrod et al.l I2OO6I ) . Ices are probably 
only present in the cores. We set out to test the desorption 
efficiencies, by including in an existing model of dark cloud 
chemistry three desorption mechanisms (desorption result- 
ing from H2 formation on grains, direct cosmic ray heating 
and cosmic ray photodesorption) that have been formulated 
in the literature, to find their relative importance in dark 
cores. We believe these three mechanisms are likely to be 
the most important in this situation, in the absence of any 
nearby X-ray sources in the molecular cloud. 

We note that there already exist several other studies 
wh ich model thes e three desorption mechanisms, for exam- 
ple IWillacv et ail l|l994bl ) (hereafter WRW94) performed a 
study very close to our own. However, our work differs to 
that of WRW94 because we try to take a simpler approach; 
by modelling a static cloud rather than a collapsing one we 
can easily see the effects of adding desorption to our chem- 
ical model. We also look further into the assumptions made 
about the ability of each desorption mechanism to desorb 
molecules with higher adsorption energies, such as H2O and 
NH3 (Section I4.3[) . Most importantly, we attempt to con- 
strain the efficiency of each mechanism by comparing our 
results to observations of CO depletion in star-forming re- 
gions. 

The paper is structured as follows: In Section [2] we give 
a brief summary of these desorption mechanisms and in Sec- 
tion [3] we describe the model. The results are given in Sec- 
tion [4] and are discussed by comparing them to existing 
models in Section [S] and observations in Section HJ] Conclud- 
ing remarks are given in Section [7] 



2 DESORPTION MECHANISMS 

2.1 Desorption resulting from H2 formation 

It has long been suggested that the energy released from 
exothermic reactions on grain surf aces can release energy ca - 
pable of desorbing mantle species (|Allen fc R obinson 19751 ). 
In particular, laboratory experiments on graphite substrates 
suggest that a non-negligible f raction, perhaps up to 40% 
l|Creighan. Perry fc Price] 120061 ). (F. Islam, private commu- 
nication) of the ~ 4.5 eV released in the surface formation of 
molecular hydrogen is transfered to the grain surface. This 
leads to local heating which may thermal l y des orb weakly 
bound mantle species (|Dulev fc Williams! Il993l ). although 
the extent to which the temperature rises transiently is not 
well determined, and depends on the local conductivity. In 
amorphous materials, the temperature rise may be relatively 
large. 

Such a process may be selective, depending on the tem- 



perature ac hieved as a result of the energy deposition. In p re- 
vious work l|Dulev fc Williamslll993l : IWillacv et al ]|l994al lbl), 
the conservative assumption was made that only the most 
volatile species (such as CO, N2, NO, O2, C2 and CH4) 
would be desorbed during the transient heating, correspond- 
ing to a threshold adsorption energy of E t = 1210 K. Ini- 
tially, we shall adopt the same threshold and range of volatile 
species, but we shall also consider the effects of variations 
from that value in Section |4, 31 

The rate of desorption by this process for species i is 
given by: 



Rm = sRu 2 M B (i,t) cm 3 s 



(1) 



where Ru 2 is the rate of H2 formation on grains. Rn 2 is pro- 
portional to the available grain surface area and the stick- 
ing/reaction probability to form H2 and is empirically con- 
strained (at T=10 K) to be given by: 



i?H 2 = 3.16 x 10~ 17 n(H)n H cm" 3 s" 



(2) 



where n(H) is the number density of atomic hydrogen and 
7ih is the total hydrogen nucleon density. M s (i,t) is the frac- 
tion of the mantle consisting of species i, calculated self- 
consistently as a function of time, e is an efficiency param- 
eter such that eM s (i,t) gives the number of molecules of 
species i desorbing per H2 molecule formed. The efficiency 
of this process is uncertain, but we can make a rough up- 
per estimate as follows: If we assume that 40% (~ 2 eV) of 
the energy released during each H2 formation is transfered 
to the grain, and given that the adsorption energies of the 
species we consider are ~ 1000 K (~ 8.2 x 10~ 2 eV), wo 
would expect that a maximum of ~ 20 molecules could be 
desorbed every time an H2 molecule is formed. It is likely, 
however, that this process is less efficient, because such a 
high value for e would prevent the build-up of any m antle 
material until very late times l|Willacv fc Millarlll998 l ) . We 
therefore run our model with values of e ranging from 0.01 
to 1.0, and present results for e = 0.01 and 0.1. 

2.2 Desorption by direct cosmic ray heating 

In this paper, the rate of desorption by direct cosmic ray 
heating is calculated by simply considering the number of 
molecules capable of being desorbed per cosmic ray impact. 
This rate is different to that adopted in other models (see 
the discussion in Section r5.4|l . 

As with desorption via H2 formation, th is process is 
believed to be selective (I Willacy et al.|[l994bl ). so only the 
volatile species (CO, etc.) are expected to be desorbed, at a 
rate given by: 



-Rcr = F cr {iraln g )</>M s (i,t) cm 



(3) 



where _F cr is the flux of cosmic rays, <f> is an efficiency param- 
eter such that <j)M s (i,t) is the number of molecules released 
per cosmic ray impact, a g is the grain radius and n g is the 
number density of grains. We have used a value for the total 
grain surface ar ea per cm , (naln K ), o f 2.4 x 10~ 22 hh cm , 
consistent with iRawlings et alj (|l992l ) if we use a value for 
the depletion co e fficien t of ~ 1.0. 

iLeger et all (|l985l ) showed that the iron nuclei compo- 
nent of cosmic rays should be the most effective in heat- 
ing the grain, and therefore we adopt for the value for 
F CT of 2.06 x 10 -3 cm -2 s _1 which is an estimate of the 
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iron fraction of the canonical cosmic ray flux. It is based 
on the p roton cosmic ray fl ux in the interstellar medium 
given by Mor fill et al.l l|l976 ), and the i r on to proton ratio 

of 1.6 x 10" 4 (estimated by linger etail (Il985hj 

An estimate for cj> can be taken from lLeeer et all (|l985l ) 
who worked out the evaporation rate of CO per grain result- 
ing from both whole-grain and spot-grain heating by cosmic 
rays. For spot-grain heating, they calculated that 6 x 10 4 
CO molecules would be released per cosmic ray impact, in- 
dependent of grain size. For whole grain heating, however, 
the result was found to strongly depend on grain size, dom- 
inating spot-grain heating for a g ^ 2500 A. We have tested 
this mechanism with values of 4> ranging from 10 2 to 10 6 , but 
we only present the results for cj> — 10 to be consistent with 
the lower estimate predicted by spot-grain heating. Since 
the volatile species make up approximately 40% of the grain 
mantle, using <f> = 10 implies that the number of molecules 
released per cosmic ray impact is 4>M B (i) ~4x 10 4 . 



2.3 Cosmic ray induced photodesorption 

As cosmic rays travel throug h a molecular core, they io nise 
and excite the absorbing gas. IPrasad fc Tarafdarl l| 19831 ) in- 
vestigated the excitation of the Lyman and Werner systems 
of the hydrogen molecule, caused by collisions with either 
primary cosmic ray particles or secondary electrons released 
by cosmic ray ionisation. They found that emissions result- 
ing from these excitations may be capable of maintaining a 
chemically significant level of UV photon flux in the inte- 
riors of dark clouds, where the interstellar UV radiation is 
heavily extinguished (Av ^ 5). 

When a UV photon impinges upon the mantle of a 
grain, it is most likely to be absorb ed by H2O which may 
then be dissociated into H and OH (|Hartquist fc Williams! 
1990). The energy of the dissociation products cause local 
heating, capable of desorbing nearby molecules. This process 
is non-selective, and proceeds at a rate given by: 



i?cr P d = F P (naln g )YM s (i,t) 



(4) 



wher e Fp is the photon flux a nd Y is the yield per pho- 
ton. IPrasad fc Tarafdarl (|l983l ) estimated fp to be just 



but here we use a value of 4875 cm 



1350 cm"" 

taken from equation (21) of ICecchi-Pestellini fc Aiellol 
il 19921) F|. This value is m uch higher than the estimate of 
IPrasad fc Tarafdarl (|l983t ) because it explicitly includes the 
contribution to the excitation rate of H2 by cosmic ray pro- 
tons, not just the secondary el ectrons. We present the result s 
for Y — 0.1, as estimated by lHartauist fc William si |l990), 
but as this value is uncertain we also ran the model with 
values of Y ranging from 1 x 10~ 3 to 100. 

Obviously, the relative importance of these mechanisms 
depends upon the values used for e, <j> and Y, which can 
be treated as free parameters. In this paper, we present the 
results for the values stated above, and in Section[6]we derive 
upper limits for these parameters based on observations. 



1 This is consistent with using a cosmic ray ionisation rate, f , of 
1.3 X 1CT 17 s" 1 . 



3 THE MODEL 

In order to see clearly the relative importance of the vari- 
ous desorption mechanisms for the interstellar chemistry of 
molecular cores, we have used a simple one-point model of a 
static, isothermal dark cloud, at a density of rh = 10 5 cm -3 
and a temperature of 10 K. The core has a visual extinction, 
Av, of 10 magnitudes, and a radius of ~ 0.05 pc, implying 
a mass of 1 M©. 



3.1 The chemical model 



The model, adapted from IViti et all ([2004) , includes gas- 
phase reactions, freeze-out, surface reactions and desorption, 
with 127 gas phase and 40 mantle species. The gas phase re- 
action network co nsists of 1741 reactio ns from the UMIST 
Rate99 database \he Teuff et all 12000). Photoreactions are 
included, which take into account both the external inter- 
stellar radiation field and the internal cosmic ray induced 
UV field. Both direct and indirect ionisation by cosmic rays 
are also included, using a cosmic ray ionisation rate, £, of 
1.3 x 10 -17 s _1 . The hydrogen atom abundance (that drives 
the H2 formation desorption mechanism) is determined self- 
consistently in the chemical network. T he species and re - 
actions are taken from the model as in IViti et al.l l|2004l) . 
which describes the rich chemistry of hot cores. The num- 
ber of species and elements should therefore be more than 
sufficient to describe the low mass case in this paper. Freeze- 
out and grain chemistry are described in Section \3. 2 1 

The initial conditions are atomic, apart from carbon 
which is all singly ionised, and half the hydrogen nuclei are in 
the f orm of H2. The initia l elemental abundances are taken 
from lSofia fc Meveri (l200ll ). 



3.2 Freeze-out and grain chemistry 

The rate of accretion onto grains of species i (in cm~ 3 s _1 ) 



is given by IjRawlings et al 



grams 
1992): 



dn(i) 
dt 



-4.57 x 10 4 (■Raln e )T 1/2 CS(i)m(i)- 1/2 n(i) (5) 



where T is the gas temperature, m(i) is the mass of species 
i in atomic mass units, n(i) is the number density of species 
i, (irag-rig) has the same definition as in Section [2] and S(i) 
is the sticking coefficient of species (£). Following WRW94, 
we have used S(i) = 0.3 for all species. 

C is a factor which takes into account electrostatic ef- 
fects: 

C =1 for neutral species 

= 1 + (16.71 x 10 _4 /a g T) for singly charged (6) 

positive ions. 

The only surface reactions that occur are hydrogenation 
of certain unsaturated species, and dissociative recombina- 
tion of molecular ions. These reactions are assumed to occur 
instantly and the products remain on the grain surface un- 
til they are desorbed by one of the processes named above. 
H2 is an exception in that it is assumed to desorb immedi- 
ately because of the high exotherm icity of the reaction, as is 
indic ated by experimental studies (|Creighan. Perry fc Price! 
l2006h . 
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4 RESULTS 

By comparing the rate coefficients given for each of the des- 
orption mechanisms, if we use e = 0.01, <f) = 10 and Y = 0.1 
at a density of nu — 10 5 cm -3 , desorption via H2 formation 
will dominate if the equilibrium atomic hydrogen density is 
greater than ~ 0.4 cm -3 . If e is as high as 0.1 then this 
mechanism will always be dominant. 

Figures [1] and [2] compare the effect of the different 
desorption mechanisms in a static cloud at a density of 
10 5 cm -3 for selected observable molecules. In general, des- 
orption has the most significant effect on gas phase abun- 
dances after approximately 10 6 years, which is when the 
molecules begin to freeze-out. Each of the desorption mech- 
anisms are effective enough to be able to compete with 
freeze-out at late times, preventing full depletion of species 
onto grains. The exceptions are the sulphur bearing species, 
H2S and CS, whose freeze-out is only prevented by the non- 
selective cosmic ray photodesorption. This mechanism is also 
able to enhance abundances of molecules such as NH3 and 
H2S at times as early as 10 years. For the case of NH3 
this is because the NH3 in the mantle builds up relatively 
quickly (compared to molecules like CO), so direct desorp- 
tion by this mechanism can proceed close to its maximum 
rate. However, since the H2 formation and direct cosmic ray 
heating mechanisms are not able to desorb NH3 directly, 
desorption via these mechanisms can only affect NH3 once 
enough NO has been desorbed to enhance the gas phase 
production (via the reaction NO + NH+ -> NH 3 + NO+). 

For most of the species shown in Figures [T] and [2j des- 
orption via H2 formation with e = 0.1 is the most effective 
mechanism, being able to produce equilibrium abundances 
up to an order of magnitude greater than the other mecha- 
nisms. However, there are several molecules (NH3, H2S and 
CS) , for which cosmic ray induced photodesorption is more 
effective, even though desorption via this mechanism pro- 
ceeds at a rate more than ten times slower than H2 forma- 
tion with e = 0.1. This is because we have assumed desorp- 
tion via H2 formation is selective, only being able to desorb 
molecules with adsorption energies less than 1210 K (CO, 
NO, N2, O2, C2 and CH4). Evidently, the higher gas phase 
abundances of these molecules caused by selective desorp- 
tion do not have a significant effect on the gas phase chem- 
istry of molecules like H2S and CS, so their freeze-out is not 
prevented by selective mechanisms. 

In the following subsections we look at the effectiveness 
of each mechanism under different conditions, by varying 
the density and the initial atomic hydrogen density. We also 
discuss further the effects of selectivity in Section \4. 31 



4.1 Varying the density 

When the model was run at a density of riH = 10 6 cm~ 3 , it 
was found that the equilibrium abundances of the desorbed 
species were reduced approximately by a factor of 10 com- 
pared to the equivalent runs at 10 5 cm -3 . Given that n(H) 
and M s appear to show little variation in their equilibrium 
values at different densities, the rates for each desorption 
mechanism should not depend on density. However, the rate 
for freeze-out has a direct linear dependence on density, so 
it is expected that the equilibrium abundances of desorbed 



species should scale roughly as n H , which would explain the 
reduction in equilibrium abundances found in the results. 

4.2 Varying the initial atomic hydrogen density 

Since desorption via H2 formation depends on the abun- 
dance of atomic hydrogen, it would be expected that vary- 
ing the initial ratio of atomic to molecular hydrogen could 
also influence how efficient this mechanism is. However, even 
when the model was run with all hydrogen in the form of 
H2 initially, after 10 years the abundance of atomic hydro- 
gen had reached the same value (n(H) ~ 0.4 cm -3 ) as in 
the previous runs (which had only half the hydrogen nuclei 
in H2 initially), so there was no noticeable difference in the 
desorption via H2 formation rate. 

4.3 Varying the threshold adsorption energy for 
selective mechanisms 

So far we have assumed that the selective desorption mecha- 
nisms (H2 formation and direct cosmic ray heating) are only 
capable of desorbing molecules with adsorption energies less 
than or equal to the value E t = 1210 K. This number was 
chosen to be consistent with the study by WRW94, who 
assumed that only CO, N2, C2, O2 and NO could be des- 
orbed by these processes. CH4 also has an adsorption en- 
ergy of less than 1210 K, so we have also inclu ded CH4 as 
a 'volatile' species. Although iLeger et al.l (|l985l ) calculated 
that mantles composed purely of refractory ices such as CO2, 
H 2 CO, HCN, NH 3 and H 2 are unlikely to be affected by 
these mechanisms, they predicted that if these molecules are 
mixed with volatile species then spot heating processes may 
be able to raise the temperature of the grains enough to des- 
orb these refractory molecules. Since the H2 formation and 
direct cosmic ray heating mechanisms are both capable of 
spot-heating, we thought it was necessary to investigate the 
effect of varying E t . 

Figures [3] and [4] show the effect of varying E t for the H2 
formation mechanism with e = 0.01. The ad sorption ener- 
gies fo r each mantle species were taken from lAikawa et al. 
(1997), apart from H2O which was taken from lFraser et al. 
(2001). Some of the molecules, such as CO, CH4 and H2CO 
are relatively unaffected, even when E t is so high (10 000 K) 
that all molecules can be desorbed. In the case of H2CO 
which has an adsorption energy of 1760 K ( Aika wa et al.l 
Il997h this may seem rather surprising. The behaviour can 
be explained by the fact that, for the model and parame- 
ters that we have considered, gas-phase formation of H2CO 
(via CH3+O) dominates over desorption. On the other hand, 
molecules such as NH3 and HCN show differences in abun- 
dances of more than one order of magnitude when they reach 
equilibrium. These differences can be explained as follows: 

(i) NH3: The equilibrium fractional abundance of NH3 
increases from 1.4 x 10 -9 to 2.5 x 10 -8 when E t increases 
from 3000 K to 4000 K. This change is simply because the 
adsorption energy of NH 3 is 3080 K, so for E t > 4000 K, 
NH3 can be desorbed directly. For E t ^ 4000 K, the main 
production of NH3 is through the reaction NH^ + NO — » 
NO + + NH 3 . 

(ii) HCN: The equilibrium fractional abundance of HCN 
increases from 9.5 x 10~ n to 2.0 x 10 -9 when Et increases 
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Figure 1. The time evolution of the molecules NH3, H2S, CH4, CO and O2 in a static cloud of density rtpj = 10 cm" , T = 10 K, 
at a point with Ay = 10 mag. The abundances of species relative to hydrogen are shown as a function of time in years. The different 
curves compare the evolution of the species with no desorption, desorption via H2 formation with e = 0.1 and 0.01, desorption via direct 
cosmic ray heating and cosmic ray photodesorption (see key). 



from 3000 K to 4000 K. The adsorption energy of HCN is 
actually greater than 4000 K, so this is not a consequence of 
direct HCN desorption. In fact, it is the direct desorption of 
NH 3 that produces the HCN via the reaction NH 3 + CN — > 
HCN + NH 2 when E t is greater than 4000 K. Below this 
energy, the main route of production of HCN is the reaction 
H + H 2 CN^HCN + H 2 . 

For CS and H 2 S, freeze-out is only prevented if E t is greater 
than 2000 K. This is because at this energy both CS and H 2 S 
can be directly desorbed. For CS if Et is increased further 
from 2000 K to 3000 K, the equilibrium abundance of CS 
increases again by approximately an order of magnitude. 
This is due to the direct desorption of H 2 CS at 2250 K, 
which can increase the CS abundance through the reaction 
C+ +H 2 CS -» CS + CH+. 



5 COMPARISON WITH EXISTING MODELS 
AND FORMULATIONS 

We first of all compare our results to those of WRW94, 
which is the closest previous study to our own, then we 



consider two other studies, IWillacv fc Millarl (l998) and 
iHartauist fc Williams! <|l990l ) (hereafter WM98 and HW90 
respectively), whose conclusions differed to ours. We also 
discuss various assumptions that are made in most chemical 
models of molecular clouds in which desorption is included, 
with particular reference to the formulation of desorption by 
direct cosmic ray heating by HH93. 



5.1 Comparison to WRW94 

WRW94 investigated each of the desorption mechanisms dis- 
cussed above, in the case of a collapsing cloud. They con- 
cluded that the only significant mechanism was desorption 
arising from H 2 formation on grains if it was assumed to 
be non-selective, using values for e of 0.1 and 0.8. Since we 
have found that, for a static cloud, all desorption mecha- 
nisms are significant and can prevent complete freeze-out, it 
is therefore worth checking that our model agrees with that 
of WRW94 when it is run under similar conditions. 

To simulate a collapse, as in WRW94, we have used a 
spherically symmetric, isothermal model, undergoing a free- 
fall collapse. The collapse is modified by a retardation factor 
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no desorptio'n 
H2 formation:epsilon = 0. 1 
H2 formation:epsilon = 0.01 
CR heating 
CR photodesorption 




Figure 2. As in Figure [T] but for the molecules HCO+, H 2 CO, HCN, CS and N 2 H+. 



B du e to magnetic fields and other factors (Rawlings ct al. 
119921 ). WRW94 used B = 0.7, a temperature of 10 K, an 
initial density, no, of 3 x 10 3 cm -3 and a final density of 10 6 
cm -3 . The initial visual extinction of the cloud was 4.4 mag- 
nitudes. Figure [5] shows the results from our model when it 
was run with these parameters. We have plotted the abun- 
dances of the same molecules that appear in Figure 1 of 
WRW94, and have converted time to radius as they did to 
allow for easy comparison between our results and theirs. 
Even though the chemical reaction network in our model is 
much more complicated than in the WRW94 model, both 
models agree that desorption has the greatest effect on NH3 
and HCN. Also, under these conditions, we agree that the 
other mechanisms (direct cosmic ray heating and cosmic ray 
photodesorption) do not significantly affect the abundances 
of these molecules. 



5.2 Comparison to WM98 

WM98 carried out an investigation in which they tested the 
three desorption mechanisms used in our model, as well des- 
orption by the explosive chemical reactions of radicals. Their 
model used similar physical conditions to ours, except that 



they adopted a lower density (hh = 2 x 10 cm 



They 



found that models that include desorption by H2 formation 



(with values of e ranging from 0.1 to 1.0) gave the best agree- 
ment to the observed gas phase abundances in TMC-1. 

This actually agrees well with our model, since we also 
find that desorption by H2 formation with e = 0.1 is the 
most efficient mechanism. WM98 concluded, however, that 
desorption by direct cosmic ray heating and cosmic ray 
photodesorption were too inefficient because they generally 
under-predicted the gas phase abundances. Thus, for exam- 
ple, CO, was found to be underabundant by approximately 
two orders of magnitude at an age of 5 Myr using the cosmic 
ray induced mechanisms. We agree that, when comparing 
the desorption mechanisms using our standard values of e, 
<j> and Y, desorption by H2 formation can indeed maintain 
gas-phase abundances up to two orders of magnitude higher 
than the other mechanisms. But, since these parameters are 
so uncertain, we believe it is too early to confidently rule 
out cosmic ray heating and cosmic ray photodesorption as 
important desorption mechanisms. 



5.3 Comparison to HW90 

In their study, iHartauist fc Williams! (|l990l ) carried out a 
purely theoretical analysis of cosmic ray induced desorp- 
tion. They argued that direct cosmic ray heating would not 
be able to maintain significant gas-phase CO abundances be- 
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Figure 3. The effect of varying Et for desorption via H2 formation with e 
K, for the molecules NH3, H 2 S, CH 4 , CO and O2 



0.01, in a static cloud of density njj = 10 cm at 10 



cause the desorbed CO would be destroyed by reactions with 
He + (leading to the formation of hydrocarbons). Their an- 
alytical calculation of the gas-phase CO abundance at high 
densities gives the result 

n 9 (CO) oc e- 0t , 



where (3 depends on the rates of the reactions 

He + + CO -> C + + O + He, 

CO(g) ^CO(s), 

CO(s) ->CO(g) 1 

He + CR -> He+ + e + CR 

and 

He + + H 2 -> H + + H + He. 



(7) 

(8) 
(9) 
(10) 

(11) 
(12) 



Reactions © to (|10|) are the main formation and de- 
struction routes for CO, with reactions ((Sj and (|10[1 referring 
to the freeze-out and desorption of CO. Reactions (|11[) and 
(|12p are assumed to be the main formation and destruction 
rout es of He + at high dens i ties. 

lHartquist fc Williamsl (Il990l ) calculated /3 to be 



However, when we redo the calculation with the (up- 
dated) rate coefficients for these reactions used in our model, 
we obtain values of that are a factor almost six orders of 
magnitude larger. This implies that significant gas-phase CO 
abundances can be maintained over the lifetime of a typical 
cloud. The difference arises because we have used a much 
smaller desorption rate for CO (given by equation (O). So, 
in fact, using a smaller desorption rate allows for a more 
controlled release of CO into the gas phase and thus enables 
the CO to remain in the gas-phase for a longer period of 
time^ 

lHartquist fc Williamsl (|l990h also concluded that cos- 
mic ray induced photodesorption would only be effective in 
regions where the number density is 10 3 cm -3 or lower. They 
calculated the ratio at which heavy molecules desorb by cos- 
mic ray induced photodesorption to the rate at which they 
freeze-out to be 



R : 



10 5 cm- 



c 



10- 



(14) 



3 x 10 6 yr(riH/10 6 cm" 3 ). 



(13) 



where Xz is the fractional abundance of heavy molecules in 
the gas phase and £ is the cosmic ray ionisation rate. This 
expression implies that cosmic ray photodesorption will be 
more efficient at low densities. The rates we have used are 
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Figure 4. As in figure[3]but for the molecules HCO+, H 2 CO, HCN, CS and N 2 H+ 



consistent with this ratio, but our results clearly indicate 
that the cosmic ray induced photodesorption mechanism 
can still maintain significant gas phase abundances of heavy 
molecules at densities as high as 10 5 cm -3 , even though the 
efficiency of the process is not at its peak at this density. 



5.4 Discussion of other formulations and 
assumptions 

As mentioned earlier, other previous models that have in- 
cluded desorption in dark clouds tend to only include des- 
orption by direct cosmic ray heating. One of the most com- 
monly used formulations is that of HH93, who assumed that 
classical (~ 0.1 (im) grains are impulsively heated by rela- 
tivistic nuclei with energies 20-70 MeV, depositing an av- 
erage of 0.4 MeV per impact, raising the local (hot spot) 
temperature of the grain to 70 K. The subsequent (thermal) 
desorption rate is then calculated from the binding energies 
of the various molecular species and the cooling profile of 
the dust grain. In practice, the details of the cooling are 
simplified to a 'duty cycle' approach, so the desorption rate 
is proportional to the fraction of time spent by the grains 
in the vicinity of 70 K, /(70 K) ~ 3.16 x 10~ 19 . The rate 
coefficient is then given by: 



= /(70 K)n>(i) exp[-£ a (i)/70 K] s" 1 (15) 



where vo(i) and _E a (i) are the characteristic adsorbate vi- 
brational frequency and adsorption energy for species i re- 

spe ctively. 

iFlower et all (2005) formulated a similar (but simpler) 
rate for cosmic ray desorption (in cm~ 3 s~ 1 ) given by: 



R 



■crd-FPWO'i 



M g (i) (ira s n g )y exp 



-(E a (€) - Sa(CO)) 



70 K 



(16) 



where 7 = 70 cm 2 s 1 gives the desor ption rate of CO pe r 
unit area of dust grains as derived by iLeeer et all l|l985h . 
This factor takes into account the cosmic ray flux, F cr , and 
the efficiency parameter, <j>, which were explicitly included 
in our rate. The above rate includes the same exponential 
factor as in HH93, but does not take into account variations 
in fo for each species. 

We decided not to use the above formulations because, 
bearing in mind the exponential sensitivity of the desorption 
rate to the ratio of the hot spot temperature to the bind- 
ing energy, the resulting rates are extremely uncertain and 
strongly dependent on a number of very poorly constrained 
free parameters including the grain size and morphology, 
the molecular binding energies, the rate and energy/mass 
spectrum of the cosmic rays. 
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Figure 5. The variation of NH3, CH, HCO + and HCN with radius for a cloud c ollapsing from n H = 3 X 10 3 cm 3 to 10 6 cm 3 at 10 K 
including non-selective desorption resulting from H2 formation, as in Figure 1 of I Willacy et al.l dl994bl) 



We can examine a few of these assumptions in a little 
more detail: 

(i) The formulations given above assume that interstellar 
dust grains can be thought of as homogeneous, symmetri- 
cal (spherical, spheroidal, or cylindrical) entities. In reality, 
we know that the grains may have very complex morpholo- 
gies ('fluffy aggregates') with non-uniform thermal proper- 
ties and composition. We can speculate that these grains 
therefore consist of poorly thermally connected sub-units, 
so that both the spot-heating and the whole grain heating 
rates may be quite different to what is derived for spherical, 
uniform dust grains. 

(ii) Most of these studies also only consider the gas-grain 
interactions with large (~ 0.1 fiia), classical grains. The jus- 
tification given for this is that the equilibrium temperature 
of the (cosmic ray heated) smaller grains is too high to allow 
gas-phase species to freeze-out . Howe ver, this differentiation 
is not apparent in lLeger et al . (1985) and there is no direct 
observational evidence of the universal presence of a popula- 
tion of warm dust grains. The smaller grains present a much 
larger surface area per unit volume of gas than the classical 
grains, so this is an important issue. 

(iii) Concerning the desorption process itself, the assump- 



tion is usually made that the mantle ice is pure, so that the 
adopted binding/desorption energies correspond to the pure 
substance. Real ices are likely to have complex compositions 
and temperature-dependen t morphologies. Laboratory work 
l)Collings et al J 12003 . l2004h shows that the desorption char- 
acteristics are drastically modified in such circumstances in 
a way that is very species-dependent. Thus, for example, 
in a mixed CO/H2O ice, CO is desorbed in four distinct 
temperature bands (with T ~ 25 — 100 K). 

(iv) iBringa fc Johnsonl l|2004 ) performed detailed molec- 
ular dynamics calculations of the thermal evaporation pro- 
cess driven by a cylindrical heat pulse, concentrating on a 
classical grain of ~ 0.1 /jm, heated by a heavy cosmic ray 
ion. Even with the assumption of a single desorption en- 
ergy/threshold for each species they still obtained very dif- 
ferent results to HH93 - most notably the derived H2O and 
CO desorption rates were found to be an order of magni- 
tude higher and lower, respectively, as compared to HH93. 
However, even here, the calculations make very general as- 
sumptions about the cosmic ray energy spectrum and adopt 
a spherically symmetric grain with uniform thermal proper- 
ties. 

(v) The issue of the surface chemistry is often not consid- 
ered in many of the chemical models which use these des- 
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orption formulations. This would only be valid if one makes 
the seemingly unreasonable assumption that the surface res- 
idence timescale of adsorbed species is less than the surface 
migration timescale for reactive species such as hydrogen 
atoms. 

Because of these various uncertainties we must conclude 
that although approximate quantifications of the rate of des- 
orption driven by cosmic ray grain heating are useful, the 
numerical details are very poorly constrained and anything 
more complicated or specific than our simple approach is 
very hard to justify. 

The most important conclusion of this section, there- 
fore, is that due to the considerable uncertainties in the mi- 
crophysics of ice desorption which are, as yet, very poorly 
constrained by laboratory and/or theoretical studies, it is 
very difficult to make qualitative distinctions between differ- 
ent desorption mechanisms: Each of the desorption models 
discussed above has its own merits, but that there is insuf- 
ficient information to discriminate or validate these mod- 
els. However, we also make the point that, through empir- 
ically deduced values of the depletions, we can make well- 
constrained quantitative estimates of the desorption efficien- 
cies. In the next section, we compare our models to obser- 
vational results. Thus, by consideration of the deduced des- 
orption efficiencies we are able to comment on the allowed 
values of the free parameter in each of the desorption models 
(e, (f> or Y) that are consistent with the observations. 

6 COMPARISON WITH OBSERVATIONS 

As a consequence of our poor understanding of the theory of 
ice desorption, it is proposed that, rather than using desorp- 
tion efficiencies as an input to the chemical models we invert 
the process and use whatever observations we can safely in- 
terpret to constrain, empirically, the nature and efficiencies 
of the desorption processes. We can then use that infor- 
mation to predict the chemical behaviours of other species 
and generate physical diagnostic indicators of the molecular 
clouds. 

Observations of CO isotopomers in several pre-stellar 
cores provide strong evidence for the depl etion of CO onto 
grains. For example, I Redman et al.l (|2002l ) were able to es- 
timate that at least 90 % of CO is depleted within 5000 AU 
of the centre of the pre-stellar core L1689B, by comparing 
their observations of the C ir O J = 2 — ■> 1 line with mod- 
els which included CO depletion in the centre of the core. 
Since L1689B is a typical pr e-stellar core with a d ensity of 
nn = 1.2 - 1.4 x 10 5 cm" 3 (|Bacmann et al J 1200(f ). it is an 
ideal candidate with which we can compare our results. 

CO depletion has also been estimated in several other 
pre-stellar co res, such as L1544, L 1709A, L310, L328, L429 
and Oph D (|Bacmann et al.| [2002). In these cases, the CO 
depletion is estimated by comparing the observed ratio, 
X, of the C lr O and H2 col umn densities, t o the 'canoni- 
cal' abundance determined bv lFrerking et all (l982) towards 
dark cores, X can = 4.8 x 10~ 8 . The CO depletion factor, /, 
is defined by X ca , n /X, and was found to vary from 4.5 (in 
L1689B, whic h appears to be slight ly less than the depletion 
estimated by iRedman et all (|2002T to 15.5 in L429. Since 
Acan is the value of X obtained in undepleted conditions, 
and assuming that there is no selective depletion of different 



isotopomers of CO (so the fraction of C 17 depleted onto 
grains is the same as the fraction of total CO depleted), 
then the fraction of CO depleted onto grains is given by 
1 — 1//. This implies that CO depletion ranges from 78% 
to 94% in these cores. The densities of these cores are all of 
order n(H2) = 10 5 cm -3 so it is also worth comparing these 
depletion fractions to our model. 

Although our results have indicated that including des- 
orption inhibits full freeze-out allowing abundances to reach 
an equilibrium, the percentage of CO frozen-out onto grains 
is actually greater than 98% (/ = 50) for all of the mecha- 
nisms investigated, which is consistent with the observations 
above. 

If we assume that L1689B is in a state of equilibrium, 
the values of e, <j> and Y needed to give the 90% freeze- 
out estimated in L1689B are approximately 0.5, 1.3 x 10 7 
and 5.5 respectively. We calculated these values by assuming 
that freeze-out and desorption of CO are the main processes 
governing the CO abundance, so by equating these rates 
with 90% of CO on the grains and 10% remaining in the gas 
phase we can predict the efficiencies needed. We then ran the 
model with these values of e, <j> and Y to confirm that the 
observed CO depletion is obtained. Since it is possible that 
L1689B is not in equilibrium and further freeze-out could be 
achieved, this gives us an upper limit for th ese parameters. 

If we use the values estimated by iBacmann et alJ 

(2002), for L1689B (in which they deduce that CO is de- 
pleted by 78%) we obtain even higher limits for e, <f> and Y 
of 1.4, 3.3 x 10 7 and 14.0 respectively. For L429, in which a 
CO depletion of 94% was deduced, the values for e, <f> and Y 
are reduced to 0.31, 7.5 x 10 6 and 3. 2 respectively. The val - 
ues for all the other cores studied bv lBacmann et alJ (|2002h 
would lie in between the values for L1689B and L429. 

These values are much higher than the estimates used in 
our model, indicating that if these desorption mechanisms 
operate then they are very efficient. However, using these 
observations it is not possible to determine which of the 
mechanisms is operating, or if it is a combination of all three. 



7 CONCLUSIONS 

There are two major conclusions from this study: Firstly, 
we have shown that the usual assumption that cosmic ray 
desorption is the the most effective desorption mechanism in 
dark molecular clouds is not always valid. All three desorp- 
tion mechanisms that we have considered (desorption via H2 
formation on grain surfaces, direct cosmic ray heating and 
cosmic ray photodesorption) have been shown to have signif- 
icant effects on the gas phase abundances in quiescent dark 
molecular clouds and so should not be neglected in chemi- 
cal models. These processes all operate on timescales of the 
order of ~ 10 6 years. Each of the processes is capable of 
preventing total freeze-out, but in an equilibrium quiescent 
dark cloud of density 10 s cm -3 , the predicted percentage of 
freeze-out is always greater than 98%. This figure is in good 
agreement with the observations (e.g. L1689B). 

Addressing the specifics of the desorption processes, 
desorption via H2 formation, if it is efficient (e = 0.1), is the 
most effective mechanism. However the complete freeze-out 
of some species, such as CS and H2S, can only be prevented 
by the cosmic ray photodesorption mechanism, which is non- 
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selective. The relative importance of the three mechanisms 
appear to be insensitive to variations in the density and the 
initial atomic to molecular hydrogen ratio. For the selective 
desorption mechanisms, choosing the threshold adsorption 
energy, E t , (such that molecules with adsorption energies 
less than E t will be desorbed), can have a strong effect on 
the chemistry, particularly on molecules such as NH3, HCN, 
CS and H 2 S. 

Secondly, our understandings of the chemical and phys- 
ical structure of dust grains and the physical processes 
which drive desorption are, as yet, very incomplete and a 
purely theoretical approach to the problem is inadvisable. 
In this study we have investigated three particular desorp- 
tion mechanisms (desorption via H2 formation, direct cosmic 
ray heating and cosmic ray photodesorption) and have used 
the observed molecular depletions to constrain the poorly- 
determined free parameters in desorption processes. This is 
the first attempt at an empirical determination of the des- 
orption efficiencies. Bearing in mind the huge uncertainties 
in these efficiencies, this empirical approach is the one that 
we recommend for use in future studies of interstellar chem- 
istry where gas-grain interactions play an important role. 
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